Role of backflow correlations for the non-magnetic phase of the t—t' Hubbard model 
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We introduce an efficient way to improve the accuracy of projected wave functions, widely used to 
study the two-dimensional Hubbard model. Taking the clue from the backflow contribution, whose 
relevance has been emphasized for various interacting systems on the continuum, we consider many- 
body correlations to construct a suitable approximation for the ground state at intermediate and 
strong couplings. In particular, we study the phase diagram of the frustrated t—t' Hubbard model 
on the square lattice and show that, thanks to backflow correlations, an insulating and non-magnetic 
phase can be stabilized at strong coupling and sufficiently large frustrating ratio t' /t. 

PACS numbers: 71.10.Fd, 71.27.+a, 71.30.+h, 75.10.-b 
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Introduction. Recently, the interest in the role of frus- 
trating interactions in electronic systems has consider- 
ably increased since in this regime new exotic phases may 
appear. Many experiments suggest the possibility to have 
disordered phases down to very low-temperatures (much 
smaller than what one would expect from a mean-field 
approach) or even to zero temperature. Such phases are 
generically called spin liquids. In this respect, the organic 
molecular materials k-(ET)2X, X being a monovalent an- 
ioufii^ represent an interesting example, since they show 
a particularly rich phase diagram. In the conducting lay- 
ers, ET molecules are strongly dimerized and form a two- 
dimensional (2D) triangular lattice. Since the valence of 
each ET dimer is +1, the conduction band is half filled. 
By acting with an external pressure, it is possible to vary 
the ratio between the on-site Coulomb repulsion and the 
bandwidth, driving the system through a metal-insulator 
transition (MIT). 

The minimal model to describe the physics of corre- 
lated electrons is the Hubbard model 
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where ^{ci^a-) creates (destroys) an electron with spin 

a on site i, ni^^r = c\^Ci^a^ tij is the hopping ampli- 
tude, that determines the bandwidth, and U is the on-site 
Coulomb repulsion. In this work, we focus our attention 
on the half-filled case with N electrons on TV sites and 
consider the square lattice with both nearest- and next- 
nearest-neighbor hoppings, denoted by t and t' , respec- 
tively. This model represents the prototype for frustrated 
electronic materialsji^, and, recently, it has been widely 
studied by different numerical techniques, with contra- 
dictory outcomesj^i^i^ii^ Here we present the results for 
the zero-temperature phase diagram, obtained by using 
projected wave functions, which only after considering 
backflow correlations are accurate enough to describe the 
highly-correlated regime. 

The variational approach. Variational wave functions 
for the unfrustrated Hubbard model, describing the anti- 
ferromagnetic phase, can be constructed by considering 
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the ground state \AF) of a mean- field Hamiltonian con- 
taining a band contribution and a magnetic term Haf = 
Aa_f e'*^'^^<5'J, where Sj is the x component of the 
spin operator Sj — {Sj , Sj , Sj ). In order to have the cor- 
rect spin-spin correlations at large distance, we have to 
apply a suitable long-range spin Jastrow factor, namely 
|*af) = Js\AF), with J, = exp[-^Y.^,J^^,JS^S^], 
which governs spin fluctuations orthogonal to the mag- 
netic field Aaf^ 

On the other hand, spin-liquid (i.e., disordered) states 
can be constructed by considering the ground state 
\BCS) of a BCS Hamiltonian and then applying to it 
the so-called Gutzwiller projector, \RVB) = Vg\BCS), 

where Vg = Ylii^ — S"^i,T"'j,i) ^^'^ 9 = 1 '^°'^^ In pure 
spin models, where the U is infinite and charge fluctu- 
ations are completely frozen, these kind of states can 
be remarkably accurate and provide important predic- 
tions on the stabilization of disordered spin-liquid ground 
states. ^^'^■^ However, whenever U/t is finite, the varia- 
tional state must also contain charge fluctuations. In 
this regard, the simplest generalization of the Gutzwiller 
projector with g < I, that allows doubly occupied sites, 
is known to lead to a metallic phaseJ^ In order to ob- 
tain a Mott insulator with no magnetic order, it is neces- 
sary to consider a sufficiently long-range Jastrow factor 
J = exp[-i X^i.j = J2a'^i,<^ being the lo- 

cal density. Nevertheless, the accuracy of the resulting 
wave function \'^bcs} = J\BCS) can be rather poor in 
2D for large on-site interactions,^^ especially in the pres- 
ence of frustration (see below). Therefore, other contri- 
butions beyond the Jastrow factor must be included. 

An alternative way, suitable to describe the strong- 
coupling regime, is to start from the fully-projected 
state \RVB) and then apply the unitary transforma- 
tion, that was introduced long time ago to connect the 
Heisenberg and the Hubbard models, namely s) — 
exp{iS)\RVB) . This kind of approach is rather difficult 
to implement for large clusters, since, in contrast to the 
Jastrow term, S is non-diagonal in the natural basis \x) 
where the electrons with spins quantized along z occupy 
the lattice sites. The expression valid for strong coupling. 
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can be introduced by considering a linear expansion of 
each single-particle orbital: 



FIG. 1: (Color online) Results for 18 electrons on 18 sites 
as a function of U/t. Upper panels: Accuracy of energy 
AE — {Eo — Ev), Ev and Eo being the variational and the ex- 
act values, respectively. Lower panels: Overlap between the 
exact ground state and the variational wave functions. The 
BCS state with long-range Jastrow factor is denoted by blue 
triangles, the BCS state with backflow correlations and Jas- 
trow term by red circles. The results for AE considering one 
Lanczos step upon the BCS state, i.e., (1 4- aTi.)\'i'), are also 
shown (black squares). 



i.e., (1 + iS)\RVB), is clearly not accurate for large sys- 
tem sizes, since it can allow a single doubly-occupied site 
at most. In this respect, accurate results for small clus- 
ters can be also obtained by performing one Lanczos step, 
(1 + an)\'i/BCs)r or by considering exp{hK)VG\BCS) 
(where h and g are variational parameters and K is the 
hopping Hamiltonian) 

The backflow wave function. In order to improve in 
a size-consistent way the previous wave functions |^'a_f) 
and |^'bcs)j we want to modify the single-particle or- 
bitals,^^ in the same spirit of the backflow correlations, 
which have been proposed long time ago by Feynman and 
Cohen, to obtain a quantitative description of the roton 
excitation in liquid Helium.^'' The backflow has been im- 
plemented within quantum Monte Carlo calculations to 
study bulk liquid '^Hc j^^'^^ and used to improve the de- 
scription of the electron jellium both in two and three 
dimensions. ^■^'^'^ More recently, it has been applied to 
metallic Hydrogen. ■^^ Originally, the backflow term corre- 
sponds to consider fictitious coordinates of the electrons 
r^, which depend upon the positions of the other parti- 
cles, so to create a return flow of current: 
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where are the actual electonic positions and 77a, /3[a;] are 
variational parameters depending in principle on all the 
electronic coordinates, namely on the many-body con- 
figuration \x). The variational wave function is then 
constructed by means of the orbitals calculated in the 
new positions, i.e., (/'(r^). Alternatively, the backflow 
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where CQ.^[a;] are suitable coefficients. The definition ([3]) 
is particularly useful in lattice models, where the coor- 
dinates of the particles may assume only discrete values. 
In particular, in the Hubbard model, the form of the new 
"orbitals" can be fixed by considering the U ^ t limit, 
so to favor a recombination of neighboring charge fluctu- 
ations (i.e., empty and doubly-occupied sites): 



(4) 



where we used the notation that 4>k{'^i^a) — (O|ci^o-|0fc), 
being \(j)k) the eigenstates of the mean- field Hamiltonian, 
Di = n^^fUi^i, Hi = hi^^h.^i, with h^^^ = (1 - ?ii,cr), 
so that Di and Hi are non zero only if the site i is 
doubly occupied or empty, respectively; finally e and 77 
are variational parameters (we can assume that e = 1 if 
DiHj =0). As a consequence, already the determinant 
part of the wave function includes correlation effects, due 
to the presence of the many body operator DiHj. The 
previous definition of the backflow term preserves the 
spin SU(2) symmetry. A further generalization of the 
new "orbitals" can be made, by taking all the possible 
virtual hoppings of the electrons: 

j 

'72 ^iij (rii^ahi^anj^^ahj^) (f>k{rj^a) + 

3 

^3 ^ {D.,n.j^^„h.j„ + n.,^c,hi-aHj) (j)k{vj^a), (5) 


where e, 771, 772, and 7^3 are variational parameters. The 
latter two variational parameters are particularly impor- 
tant for the metallic phase at small U/t, whereas they 
give only a slight improvement of the variational wave 
function in the insulator at strong coupling. The deflni- 
tion Eq. ^ may break the SU(2) symmetry, however, the 
optimized wave function always has a very small value of 
the total spin square, i.e., (5*^) ~ 0.001 for 50 sites. All 
the parameters of the wave function (contained in the 
mean-field Hamiltonian, in the Jastrow term, and in the 
backflow term) can be optimized by using the method of 
Ref. [1^ Finally, the variational results can be compared 
with more accurate (and still variational) ones obtained 
by Green's function Monte Carlo, implemented with 
the so-called fixed-node (FN) approximation.^^ 

Results. Let us start by considering the comparison of 
the variational results with the exact ones on the 18-site 
cluster at half filling. In Fig. [H we show the accuracy 
of the variational BCS state (with and without backflow 
correlations) and the overlap with the exact ground state 
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FIG. 2: (Color online) Variational energies per site (in unit of 
J = At^ /U) for the BCS state with long-range Jastrow factor, 
with and without backflow correlations, and 98 sites. The re- 
sults for the wave function with antiferromagnetic order and 
no BCS pairing are also shown. Arrows indicate the varia- 
tional results obtained by applying the full Gutzwiller projec- 
tion to the mean-field states for the corresponding Heisenberg 
models. 
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FIG. 3: (Color online) Lower panel: Phase diagram of the 
frustrated t—t' Hubbard model, as obtained by comparing 
the variational energies of the backflow wave functions. Upper 
panel: Comparison between the variational energies per site 
and the FN ones for U/t = l& and 98 sites. 



for two values of the frustrating ratio, i.e., t' /t — and 
0.7. The backflow term is able to highly improve the 
accuracy both for weak and strong couphngs. We also 
notice that backflow correlations are more efficient than 
applying one Lanczos step, i.e., {l+aH)\^ bcs) i that was 
used in previous calculations)^ The overlap between the 
exact ground state and the backflow state remains very 
high, even for large U, and the improvement with respect 
to the BCS state is crucial, especially in the frustrated 
regime. 

Backflow correlations remain efficient also for larger 
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FIG. 4: (Color online) Upper panel: Variational (empty sym- 
bols) and FN (full symbols) results for ^(Tr, n) divided by A''. 
Both calculations have been done by using the projected BCS 
wave function; U/t = 16 and t'/t = (triangles), U/t = 24 
and t'/t = 0.7 (circles), and U/t — S and t'/t = 0.75 (squares). 
Lines are guides to the eye. Lower panel: Variational results 
for N{q) divided by |g| for 98 (empty symbols) and 162 (full 
symbols) sites and t'/t — 0.75; from top to bottom: U/t — 4, 
6, 7, 8, and 16. 



sizes and provide much lower energy than the Lanczos 
step wave function, e.g., for 98 sites with U/t = 20 and 
t'/t ~ 0.7, the energy per site with the backflow wave 
function is Eb/t = —0.2352(1), while the one with one 
Lanczos step is Eis/t — —0.2310(1) (for 18 sites they 
are Eb/t = -0.23741 and Ei^/t = -0.23566). The FN 
energy obtained with the backflow state is Ep^/t = 
—0.2395(1), rather close to our estimation of the exact 
value (based upon an extrapolation obtained with zero 
and one Lanczos step) that is E/t ~ —0.246. 

By increasing U/t, the variational energy extrapolates 
to the one obtained by taking the fully-projected state 
\RVB) in the spin model. On the contrary, without us- 
ing backflow terms, the energy of the BCS state, even in 
presence of a fully optimized Jastrow factor, is few hun- 
dredths of J = At'^/U higher than the expected value, see 
Fig. [2] Moreover, whenever frustration is large enough, 
backflow correlations are useful also in the antiferromag- 
netic state l^'yii?), while for t' ~ they are not necessary 
to extrapolate correctly to the value of the spin model, 
see Fig. H 

In order to draw the ground-state phase diagram of 
the t—t' Hubbard model, we consider three different wave 
functions with backflow correlations: Two antiferromag- 
netic states I^I/af) with Q = (7r,7r) and Q = (tt, 0), rele- 
vant for small and large t'/t, and the non-magnetic state 
I'^BCs)- The variational phase diagram is reported in 
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Fig. [3l The first important outcome is that, without 
backflow terms, the energies of the spin-liquid wave func- 
tion are always higher than those of the magnetically or- 
dered states, for any value of frustration t' /t. Instead, by 
inserting backflow correlations, a spin-liquid phase can be 
stabilized at large enough U/t and frustration (see also 
Fig. [2). The small energy difference between the pure 
variational and the FN energies demonstrates the accu- 
racy of the backflow states, see Fig. [3] Notice that |^'a_f} 
and I'^BCs) have different nodal surfaces, implying dif- 
ferent FN energies. 

For small Coulomb repulsion and finite t' /t the static 
density-density correlations N{q) = {n-qUq) (where Uq is 
the Fourier transform of the local density rii) have a lin- 
ear behavior for \q\ — > 0, typical of a conducting phase. A 
very small superconducting parameter with d^2_y2 sym- 
metry can be stabilized, suggesting that long-range pair- 
ing correlations, if any, are tiny. By increasing U/t, a 
MIT is found and N{q) acquires a quadratic behavior in 
the insulating phase, indicating a vanishing compressibil- 
ity. This behavior does not change when considering the 
FN approach, though the metal-insulator transition may 
be slightly shifted. In Fig. 01 we show the variational 
results for N{q) as a function of U/t for t'/t = 0.75. 
The insulator just above the transition is magnetically 
ordered and the variational wave function has a large 
Aaf] the transition is likely to be first order. By further 
increasing U/t, there is a second transition to a disor- 
dered insulator. Indeed, for U/t> 14, the energy of the 
BCS wave function becomes lower than the one of the an- 
tiferromagnetic state. In this respect, the key ingredient 
to have such an insulating behavior is the presence of a 
singular Jastrow term Vq ^ 1 /q'^, that turns a BCS super- 
conductor into a Mott insulator.— In contrast to previous 



investigationsr^i^iii^ for intermediate on-site couplings, 
our calculations indicate the possibility to have a direct 
(first-order) transition between two magnetic states, see 
Fig.lSl 

In order to verify the magnetic properties obtained 
within the variational approach, we can consider the 
static spin-spin correlations S{q) — (StqSq) over the FN 
ground state. Although the FN approach may break the 
SU(2) spin symmetry, favoring a spin alignment along 
the 2 axis (this is what we find for small lattices, by a 
direct comparison with exact results), S{q) is particu- 
larly simple to evaluate within this approach;^ and it 
gives important insights into the magnetic properties of 
the ground state. In Fig. IH we report the comparison 
between the variational and the FN results by consid- 
ering the non-magnetic state I^'bcs)- Remarkably, in 
the unfrustrated case, where antiferromagnetic order is 
expected, the FN approach is able to increase spin-spin 
correlations at q — (tTjTt), even by considering the non- 
magnetic wave function to fix the nodes. A finite value 
of the magnetization is also plausible in the insulating 
region just above the metallic phase at strong frustra- 
tion (i.e., t'/t ^ 0.75), confirming the pure variational 
calculations. On the contrary, by increasing the electron 
correlation, the FN results change only slightly the vari- 
ational value of S'(7r, tt), indicating the stability of the 
disordered state. 

In conclusion, we have introduced a novel wave func- 
tion, that highly improves the accuracy of the projected 
states, used so far. Our variational ansatz is particularly 
useful to describe non-magnetic phases, that can be sta- 
bilized in the strong-coupling regime of the t—t' Hubbard 
model on the square lattice. 
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